A large X-ray flare from the Herbig Ae star V892 Tau 
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Abstract. We report the XMM- Newton observation of a large X-ray flare from the Herbig Ae star V892 Tau. The 
apparent low mass companion of V892 Tau, V892 Tau NE, is unresolved by XMM-Newton. Nevertheless there is 
compelling evidence from combined XMM-Newton and Chandra data that the origin of the flare is the Herbig 
Ae star V892 Tau. During the flare the X-ray luminosity of V892 Tau increases by a factor of ~ 15, while the 
temperature of the plasma increases from kT ~ 1.5 keV to kT ~ 8 keV. From the scaling of the flare event, based 
on hydrodynamic modeling, we conclude that a 500 G magnetic field is needed in order to confine the plasma. 
Under the assumptions that a dynamo mechanism is required to generate such a confining magnetic field and that 
surface convection is a necessary ingredient for a dynamo, our findings provide indirect evidence for the existence 
of a significant convection zone in the stellar envelope of Herbig Ae stars. 
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1. Introduction 

Understanding the genesis and early evolution of inter- 
mediate mass stars is a fundamental problem in studies 
of star formation. Owing to the differences in stellar and 
circumstellar physics, as well as in time scales, the evo- 
lution of intermediate pre-main-sequence stars is qualita- 
tively different from that of lower- and higher-mass stars. 
The past 10 years or so have witnessed an increased in- 
terest in the subject, and a key issue is the nature and 
evolution of Herbig Ae/Be stars. These pre-main-sequence 
stars, with masses ranging between about 2 and 10 Mq, 
are the more massive counterparts of T Tauri stars. They 
share with lower mass T Tauri stars peculiarities such as 
infrared excess emission, conspicuous (optical and UV) 
line emission and irregu lar photometric variability (see 

IWaters fc Waelkerislll998l for a rev iew). 

According to classical models lllbedll965j) . pre-main- 
sequence stars with masses in excess of 2 Mq are expected 
to follow fully radiative tracks once the quasi-static con- 
traction has ended. The detection of Herbig Ae/Be stars in 
X-rays ( Zinnecker fc Preibischlll994HDamiani et aJ]ll994() 
was therefore somewhat unexpected. Main sequence stars 
of the same class (from later B to later A stars) do not 
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possess either a strong stellar wind nor a corona (as no 
magnetic dynamo mechanism is available). In contrast 
Herbig Ae/Be stars appear to possess a s trong stellar wind 
llSkinner et alJ Il 99# TrWet et a,lJ ll997li and a magnetic 
dynamo associated with an outer convection zone has been 
invok ed to explain the pe riodic variation of the emission 
lines llPr aderie et alll98rj[) as well as the stellar wind itself 
l|Finkenzeller fc MundtJll984j). 

IZinnecker fc Preibischl l|l994j) proposed that the ob- 
served X-ray luminosity is linked to the stars' strong stel- 
lar winds, possibly originating in shocks due to wind insta- 
bilities and/or in the collision betw een the fast wind and 
the remnant circumstellar material. iDamiani et al 1 lll994h 
also favored a stellar wind origin for the X-ray emission 
from Herbig Be stars. Nevertheless a magnetically heated 
corona co uld also be at the origin of the observed X-ray 
emission (|Zinnecker fc PreibisclJll994 - 

The presence of winds and the possibility of magnetic 
activity make these objects interesting targets for X-ray 
studies. The unambiguous detection of flaring activity in 
any Herbig Ae/Be star would be significant as it would 
provide indirect but strong evidence for the presence of a 
magnetically confined corona and thus of an operational 
dynamo mechanism. 

So far there is only one rep orted observation of flar- 
ing activity in an Herbig Be stan lHamaguchi et alJ l|2000F ) 



performed ASCA observations or the Herbig Be star 
MCW 297, and reported the detection of a large flare 
during which the X-ray luminosity of the source increased 
by a factor of 5, with the plasma temperature increas- 
ing from 2.7 keV during the quiescent phase to 6.7 keV 
at flare maximum. The interpretation however is affected 
by some ambig uity due to the la rge ASCA point spread 
function (PSF). lTesti et all l|l99Sf) found about 20 infrared 
sources in the ASCA error circle around MWC 297, the 
majority of which are likely to be low-mass protostars. The 
peak X-ray luminosity o f ~ 5 x 10 32 erg s _1 reported by 
iHamaguchi et all l|2000|) would correspond to a very large 
- but still possible - flare for a typical l ow-mass active 
T Tauri. For example iTsuboi et alJ l|l998l) have reported 
ASCA observations of a large flare from the Weak-lined 
T Tauri star (WTTS) V773 Tau, in which the peak flare 
luminosity was at least ~ 10 33 erg s -1 . The MWC 297 
event therefore would not be exceptional for a l ow mass 
pre-main-sequence star. IHamaguchi et al . (2002b) have re- 
cently re-observed MWC 297 with Chandra, finding the 
source hundred times less luminous in X-rays than esti- 
mated from ASCA observation, and with no evidence for 
variability, thus casting some doubts on the original inter- 
pretation. 

In this article we report the observation of a large X- 
ray flare from the Herbig Ae star V892 Tau. We have 
used publicly available XMM-Newton and Chandra data 
to monitor the X-ray activity of V892 Tau and during 
one of the two XMM-Newton exposures a large flare is 
observed. While the apparent companion of V892 Tau, 
V892 Tau NE, is unresolved by the XMM-Newton PSF, 
as discussed later, the evidence that the observed flare is 
coming from the Herbig Ae star is compelling. 

The present paper is organized as follow: after a brief 
introduction of the properties of V892 Tau below, the ob- 
servations are described in Sect. [3 The spectral and tim- 
ing analysis of the data are presented in Sect. [3] and the 
results are discussed in Sec. 



1.1. The Herbig Ae star V892 Tau 

V892 Tau - also known as Elias 1 - is a young stellar ob- 
ject located in the Taurus dark cloud complex, and is sup- 
posed to be the source of the illumination for the faint neb- 
ula IC 359. The apparent m agnitude of V892 Tau is R = 
13.14 dStrom fc Stromlll 994) and its s pectral classifi cation 
; from B9 dStrom fc Stromll994]) and AO feliasl l978 | 
[ Finkenzeller fc Munddll984h to A6 jcohen fc Kuhilll979l 
iBerrilli et alJll992ilThe et al.lll994h . Esti mates for th e vi- 
sual extinct ion also vary from Ay ~ 8 l)Eliaslll978^ and 
Ay = 8.85 l|Strom fc Stromlll994l derived from a simul- 
taneous estimate of the spectral ty pe and of the apparent 
color , R- I = 1.71) to A v — 3.9 IjZinnecker fc Preibiscr] 
119941) . The star is usually placed at a distance of 140 pc 
because of its association with the Taurus dark clouds 
The bolomctric luminosity of V892 Tau is 



L ~ 38 Lq (Bcrnlli et al. 1992) and the source is variable 
in the near-infrared |Elraall978|) . 

Through nea r-i nfrare d speckle interferometry 
iKataza fc Maiharal l|l99ll ) resolved V892 Tau into 
an unresolved core and a sub-arcsec structure elongated 
in the east-west direction. They interpret the light of the 
elongated structure as being reflected within an edge-on 
circumstellar disk of moderate optical depth. Newer 
near-infrared speckle interferom etry observations were 
performed by lHaas et "all l|l997jb who favor a scenario 
in which the diffuse component is due to scattering in 
bipolar lobes with a polar axis oriented east-west. 

V892 Tau appears to be a binary system l|Leinert et alJ 

Il997h . The apparent stellar companion, hereafter referred 
to as V892 Tau N E, lies 4.1 arcsec to the northeast, at po- 
sition angle 22° 1 l|Skinner et al.ll993h . The available mea- 
surements allow a tentative classification of V892 Tau NE 
as a WTTS with spe ctral type M2 redden ed by 8-12 mag 
of visual extinction l|Leinert et alJll997|l . The study by 
iLeinert et alJ l)l997j) . based on speckle interferometry, was 
carried out with the explici t purpose of detectin g binaries 
among Herbig Ae/Be stars. ILeinert et all l)l997h achieved 
a resolution of ~ 0.1 arcsec, that at a distance of 140 pc 
correspond to 14 AU, and they do not detect other stars 
in the vicinity of V892 Tau. They conclude that in the 
near-infrared V892 Tau is basically a wide binary system. 



2. Observations 

The X-ray observations discussed in this paper were ob- 
tained with the XMM-AferatoTi and the Chandra obser- 
vatories. The XMM-Newton observations of V892 Tau 
consists of two deep (74.4 and 45.1 ks nominal) consec- 
utive exposures, the first starting on March 11 2001 at 
12:40:22 UT and the second one starting on March 12 
2001 at 10:23:10 UT. All three EPIC cameras were active 
at the time of the observations, in full-frame mode with 
the medium filters. The Principal Investigator for these 
observations is F. Walter and the observation target is 
the triple WTTS system V410 Tau. The observations are 
publicly available from the XMM-Newton archive. 

The raw XMM-Newton data have been processed by us 
with the standard SAS V5.4.1 pipeline system, concentrat- 
ing, for the spectral and timing analysis, on the EPIC-pn 
camera. In each of the two XMM-Newton exposures the 
background is affected by a large proton flares of more 
than 10 ks of duration. We have retained only time inter- 
vals in which the count rate for the whole frame of pho- 
tons above 5 keV was below a certain threshold (3.3 cts/s 
in the present case). This operation omits roughly 30% of 
the observing time, but effectively reduces the background 
level by a factor of ~ 4. Source and background photons 
were extracted using a set of scripts purposely developed 
at Palermo Observatory. 



1 In the convention adopted by Skinner et al. (1993), position 
angle 0° is north and position angle 90° is east. 



Source and background regions were denned interac- 
tively in the DS9 display software, with the background 
extracted from regions on the same CCD chip and at 
the same off-axis angle as for the source region. Response 
matrices ("rmf and ARF files") appropriate for the posi- 
tion and size of the source extraction regions were com- 
puted. The spectral analysis has been performed using the 
XSPEC package, after rebinning the source spectra to a 
minimum of 20 source counts per (variable width) spectral 
bin. 

The Chandra ACIS observation of V892 Tau was taken 
starting on March 7 2002 at 6:15:28 UT (18.0 ks nomi- 
nal). The Principal Investigator for these observations is 
P. Predehl and also in this case the observation target is 
V410 Tau 2 . The data were retrieved from the public data 
archive, with no re-processing done on the archival data. 
Source and background regions were defined in DS9, and 
light curves and spectra were extracted from the cleaned 
photon list using CIAO V. 2.2.1 threads, which were also 
used for the generation of the relative response matrices. 
Spectral analysis was performed in xspec in the same way 
as for the XMM-JVewton spectra. 

Fig-dshows the images of the V892 Tau system in the 
Chandra ACIS camera and XMM-Newton EPIC-pn cam- 
era (before and after the large flare) as well as the Palomar 
Digital Sky Survey image of the field. The 4 images are 
on the same sky coordinate scale. The source at the cen- 
tre of the images is the V892 Tau system. V892 Tau and 
V892 Tau NE are clearly resolved in the Chandra observa- 
tions. The XMM-Newton point spread function has a full 
width at half maximum of 15 arcsec and therefore cannot 
resolve the system. The source at 38 arcsec to the SE of 
the V892 Tau system is [BHS98] MHO 11, a T Tauri star 
llBriceno et all Il998l) first identified in ROSAT data by 
IStrom fc Stroml lll994h . The bright source in the North- 
East corner is Hubble 4, a well known T Tauri. 

In Tabled we report the coordinates of V892 Tau and 
V89 2 Tau NE as derived from the radio (VLA) observation 
of ISkinner et alJ l)l993|h to be compared with the coordi- 
nates of the sources in the XMM-JVew<on and Chandra 
data. The source coordinates from the EPIC-pn data cor- 
respond to the peak of a gaussian distribution fitted to the 
source image, the source coordinates for the Chandra data 
simply correspond to the brightest pixel in the source im- 
age. The agreement between the sources' positions in the 
Chandra image and as determined from VLA observations 
is excellent (within 0.4 arcsec). The source coordinates 
derived from the XMM-Newton observations before the 
flare have a 2.4 arcsec offset from V892 Tau (well within 
the uncertainty expected for the determination of posi- 
tions of EPIC X-ray sources) and a 6 arcsec offset from 
V892 Tau NE. After the flare the position offsets become 



2 There are other two Chandra exposures of V410 Tau avail- 
able in the public archive (of 15 and 11 ks respectively). In 
these observations however V892 Tau is at an extreme off-axis 
angle, and we have therefore not analysed them. 



3.1 arcsec and b.8 arcsec from V892 lau and V892 lau NE 
respectively. 



3. Spectral and timing analysis of V892 Tau and 
V892 Tau NE 

3.1. Chandra observation 

The light curves of V892 Tau and V892 Tau NE during the 
~ 18 ks Chandra observation of March 2003 are shown in 
Fig. El In the Chandra data, The Herbig Ae star V892 Tau 
is clearly resolved from the less luminous companion, the 
WTTS V892 Tau NE. 

During the Chandra observation V892 Tau shows sig- 
nificant variability, with its X-ray luminosity varying by 
a factor of 2 in less than 1 ks. The rise in the source lu- 
minosity is impulsive (the source doubles its luminosity 
in less than 1 ks) and it is followed by a slow decay. A 
proper time resolved spectral analysis of the rise phase is 
not possible with the signal to noise ratio of this data set. 
We have inspected the light curves of the source in a soft 
band (0.0-1.7 keV) and in a hard band(1.7-8 keV) and 
see evidence for a hardening of the spectrum during the 
rise phase. The decay time r ~ 1 hr is typical of X-ray 
stellar flaring events. 

The bottom panel of Fig. [21 shows the light curve of 
V892 Tau NE. The source is weak, on average 8 to 10 
times less luminous than V892 Tau. Within the uncer- 
tainty of the large error-bars the source does not ap- 
pear significantly variable. The probability of constancy 
of V892 Tau NE according to the Kolmogorov-Smirnov 
test (which measures the maximum deviation of the inte- 
gral photon arrival times from a constant source model) 
is 21%. 

The ACIS spectra of V892 Tau and V892 Tau NE are 
shown in Fig. |3| As summarised in Table |3| the spec- 
trum of V892 Tau is well fitted by an absorbed 1-T 
plasma model with N(H) = 0.83 ± 0.08 x 10 22 cm" 2 
and temperature kT = 2.10 ± 0.19 keV. The metallic- 
ity of the plasma is not well constrained by the data. 
The spectral data for V892 Tau NE have poor signal 
to noise ratio, nevertheless a fit with an absorbed 1-T 
plasma model provide a useful estimate of the absorb- 
ing column density and the plasma temperature. We de- 
rive a value of N(H) = 1.20 ± 0.18 x 10 22 cm' 2 , which 
is similar to the value derived for V892 Tau and there- 
fore consistent with the hypothesis that V892 Tau NE 
is a physical companion of V892 Tau. The derived best- 
fit value of kT = 1.04 ± 0.17 keV for the plasma tem- 
perature of V892 Tau NE is instead significantly differ- 
ent from the value derived for V892 Tau. The model de- 
pendent fluxes (as derived from the best-fit models) are 
7.1 x 10 _13 erg cm" 2 s" 1 in the band 0.55-7.50 keV, and 
5.7 x 10~ 14 erg cm" 2 s" 1 , in the band 0.67-7.50 keV, re- 
spectively for V892 Tau and V892 Tau NE. 




Fig. 1. The field centered on the system V892 Tau and V892 Tau NE in the EPIC-pn camera before the flare (top left) 
and during the flare (top right), in the ACIS camera (bottom left) and in the Digital Sky Survey (bottom right). The 
images are on the same coordinate scale. The system is resolved in the Chandra data. The source at the south-east of 
the V892 Tau system is [BH98] MHO 11, a T Tauri star, the source at North-East is Hubble 4, also a T Tauri. 



Table 1. Positions for the two components of the V892 Tau system from radio (VLA) observations (Skinner et al. 1993) 
and as derived here from XMM-Newton and Chandra images. 

Source VLA XMM" XMM ft Chandra 

RA(.T2000) RA(J2000) RA(.I2000) RA(J2000) 

Dcc(J2000) Dec(J2000) Dec(J2000) Dec(J2000) 

V892 Tau 4 18 40.60 4 18 40.67 4 18 40.67 4 18 40.63 

28 19 15.9 28 19 13.7 28 19 12.9 28 19 15.8 

V892 Tau NE 4 18 40.70 c 4 18 40.74 

28 19 19.7 c c 28 19 19.3 

a X-ray quiescent, 6 X-ray flaring, c unresolved from V892 Tau. 



3.2. XMM- Newton observations 

The light curve of the V892 Tau system derived from 
the two consecutive XMM-Newton exposures are shown 
in Fig. 0] All the gaps in the light curve except the one 
around 64 ks are due to the filtering process that we have 
applied to the raw data in order to remove the effect of 
solar proton flares. In particular roughly 10 ks at the be- 
ginning of the first exposure have been removed. The gap 
around 64 ks is due to the time difference between the 
end of the first exposure and the start of the second (one 
hour). 



During the first ~ 80 ks the light curve of the 
V892 Tau-V892 Tau NE system presents significant vari- 
ability with a time scale and amplitude similar to the one 
observed for V892 Tau in the Chandra data. During the 
last 30 ks the light curve of the system undergoes a dra- 
matic variation, in what appears to be a large flare. The 
source counts increases by a factor of ~ 8 in less than 10 
ks. The source flux then remains at around the maximum 
level until the end of the XMM-Newton exposure. 

We have performed the spectral analysis of the two 
XMM-Newton exposures separately and we summarise the 
results into two separate tables: Table for the first ex- 
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Start Time 12340 6:40: 2:914 Stop Time 12340 11:30: 2:914 

Bin time: 3600. 



6000 B000 
Time (s) 



1.2x10* 1.4x10* 



Start Time 12340 7: 5: 2:914 Stop Time 12340 11: 5: 2:914 

Fig. 2. Background-corrected light curves of V892 Tau 
(top) and V892 Tau NE (bottom) from the Chandra data, 
where the two stars are well resolved. Note the different 
vertical scale. 



posurc and Table |21 for the second exposure. The reason 
for this is that while for the first exposure two tempera- 
ture plasma models are necessary to obtain acceptable fits 
to the spectra for the second exposure one-temperature 
plasma models provide adequate fits to the spectral data. 

The spectrum of the V892 Tau system during the first 
exposure of 74 ks, corresponding to an effective integra- 
tion time of ^ 64 ks (after the proton flare filtering pro- 
cess), is shown in Fig. 03 together with the best fit ab- 
sorbed 2 temperature plasma model. The best fit values 
for the model parameters are summarised in Table El to- 
gether with the reduced \ 2 of the fit and its null hypothesis 
probability, P. A model with an absorbing column density 
of N(K) = 0.92 x 10 22 cm~ 2 , plasma temperatures kT\ = 
1.02 keV and kT 2 = 2.82 keV and a metal abundance 
Z = 0.21 provides a good fit to the integrated spectrum. 
An absorbed one temperature plasma model (as used for 
the Chandra data, which however have lower statistics) 
does not provide an acceptable description of the spec- 
trum. Nevertheless a two temperature plasma model with 
a metal abundance and two plasma temperatures frozen 
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Fig. 3. Background-corrected spectra of V892 Tau (top) 
and V892 Tau NE (bottom), from the Chandra data. The 
best-fit absorbed \-T models arc also shown. 



at the values above (Z = 0.21, kT x = 1.02 keV and kT 2 = 
2.82) provides a good fit to the Chandra data (P = 0.72) 
by letting the two emission measures and the absorbing 
column density vary - see TableEl The fitted values for the 
two emission measures are EM\ = (4.3 ± 1.6) x 10 53 cm -3 
and EM 2 = (2.2 ± 0.4) x 10 53 cm" 3 , consistent with the 
values derived from the XMM-Newton spectrum. The rel- 
ative model flux level is 7.3 x 10 _13 erg cm~ 2 s _1 (in the 
0.55-7.50 keV band). The value for the absorbing col- 
umn density derived in this way from the Chandra data, 



N(R) = (1.24 ± 0.09) x 10 22 cm" 2 is higher than the 
value derived from XMM-JVewton (and a model with the 
value derived from XMM-JVeraton does not provide an ac- 
ceptable description of the data). We note though that 
a systematically higher value for the absorbing column 
density derived from Chandra data is consistent with the 
known presence of (likely carbon-based) contamination on 
the ACIS chips. This causes additional low-energy absorp- 
tion (up to 50% near the C edge) not accounted f or in the 
current response matrices l|Plucinskv et alJl2003h . 

Given the source variability during these ~ 64 ks of 
observation, we have investigated the presence of signifi- 
cant spectral variation. We have subdivided the data into 
three intervals (the first 30 ks, the following 18 ks and 
the final 16 ks), chosen to ensure similar statistics in the 
resulting spectra. As the initial part of the observation is 



rather heavily contaminated by high background, a higher 
fraction of it was discarded, and therefore the first inter- 
val is significantly longer. The three spectra were modeled 
with an absorbed two temperature plasma model. The 
fitted values of the model parameters are summarised in 
Table|3 As it can be seen by inspecting the table, the short 
term variability observed in the system V892 Tau during 
this first XMM- Newton observation does not appear to be 
associated with significant spectral changes. 

A time resolved spectral analysis was also performed 
for the second XMM-Newton observation of V892 Tau sys- 
tem, when the large flare took place. The data were sub- 
divided into three intervals: a first 14 ks interval while 
the source is quiescent, a second 8 ks interval while the 
source counts are rising and the last 10 ks, while the 
source is at its luminosity maximum 3 . The three spec- 
tra are shown in Fig. El and their best-fit parameters are 
listed in Table |21 As explained at the beginning of this 
section an absorbed one-temperature plasma model pro- 
vides an acceptable description for all the three spectra 
derived from this exposure, so we did not attempted fits 
with two-temperature plasma models (which were neces- 
sary to obtain acceptable fits of the spectra derived from 
the first XMM- Newton exposure). In addition, the ap- 
proach to flare modelling that we present in Sect. 14.2.11 
relies on single temperature paramaterization of the X-ray 
spectrum. Nevertheless, for the spectrum derived from the 
first 14 ks interval of this exposure (corresponding to the 
quiescent phase), we verified that an absorbed two tem- 
perature plasma model with a metal abundance and two 
plasma temperatures frozen at the values derived from the 
first XMM-Newton exposure (Z = 0.21, fcTi = 1.02 keV 
and kT 2 = 2.82) indeed fits the spectrum (P = 0.08). The 
value for the absorbing column density derived in this way 
is N(U) = (0.83 ± 0.08) x 10 22 cm" 2 , somewhat higher 
than the value derived with the one-temperature plasma 
fit (JV(H) = (0.65 ± 0.06) x 10 22 cm" 2 , second line of 
Table |2J), and in better agreement with the value derived 
from the first XMM-Newton exposure. This confirms the 
lack of significant spectral variation in the V892 Tau sys- 
tem during the quiescent phase. 

As can be seen from Table during the flare the only 
significant spectral variation occurs in the plasma temper- 
ature, while the fitted values for absorbing column density 
and plasma metallicity remain essentially unchanged. The 
plasma temperature increases from kT =1.53 keV before 
the flare to kT = 8.11 during the rising phase and remains 
around that value afterwards. The model dependent flux 
density of the source, in the energy band 0.35-7.50 keV, 
goes from 6.7 x 10~ 13 erg cm~ 2 s _1 before the flare to 
1.0 x 10~ n erg cm -2 s _1 during the flare maximum. The 
peak X-ray luminosity for the flare is Lx = 2.4 x 10 31 erg 
s _1 , high but not exceptionally so for stellar X-ray flares. 



3 The segment lengths do not add up to 45 ks, because of 
the gaps in the data due to the filtering out of the solar proton 
flares. 



Bin time: 600.0 s 




Time (s) 

Start Time 11979 16: 5:48:297 Stop Time 11980 21:45:48:297 

Fig. 4. Background-subtracted light curve of the system 
V892 Tau, from the two consecutive XMM-Newton obser- 
vations. With XMM-Newton the Herbig Ae star V892 Tau 
is unresolved from the companion V892 Tau NE, never- 
theless our data indicate that this light curve is completely 
dominated by the activity of V892 Tau (see discussion). 
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Fig. 5. Background corrected spectrum of the V892 Tau 
system during the first XMM-Newton observation (before 
the flare). 



4. Discussion 

During the 18 ks Chandra observation, V892 Tau (which 
is well resolved from the apparent companion) presents 
significant variability, with its X-ray flux varying by a 
factor of 2 in less than 1 ks. This type of variability 
has never been reported before for a Herbig Ae star. 
Undoubtedly, though, it is the large variation of luminos- 
ity of the V892 Tau system (the large flare) observed with 
XMM-Newton that is most remarkable. The steep impul- 
sive rise of the source counts by a factor of 10 together 
with the impulsive rise of the plasma temperature (from 
1.5 to 8.1 keV) are all consistent with the interpretation 
of the observed variation in terms of a large coronal flare. 



lable 2. Ihe first tour lines ot the table report the best-ht spectral parameters tor the V892 lau system during the 
first XMM- Newton observation (before the flare). The parameters are derived assuming a two-temperature optically 
thin plasma model. Fits with a single temperature absorbed plasma model are unacceptable in this case. The last line 
reports the result of a spectral fit to the Chandra data on V892 Tau with a two-temperature plasma model with the 
two plasma temperatures and the metal abundance frozen at the values derived from this XMM-Newton exposure 
(first line). See text. 



Time interval 


jV(H) 


kTi 


kT 2 


EMi 


EM 2 


Z 


x 2 


P 


ks 


10 22 cm" 2 


keV 


keV 


10" 3 cm-- 1 


10 M cm _i 


Z 






0-64 (total) 


0.92 ± 0.05 


1.02 ± 0.04 


2.82 ±0.41 


5.62 ±2.61 


3.56 ± 1.21 


0.21 ±0.06 


1.07 


0.21 


0-30 


0.94±0.10 


0.98±0.10 


2.69±0.32 


2.65±2.16 


4.46±1.09 


0.51±0.20 


0.92 


0.66 


30-48 


0.89±0.07 


1.08±0.09 


3.17±1.26 


7.42±4.98 


3.03±3.08 


0.13±0.09 


1.14 


0.15 


48-64 


0.98±0.13 


0.77±0.08 


1.95±0.17 


3.79±2.84 


4.43±1.11 


0.31±0.15 


1.17 


0.16 


0-18* 


1.24±0.09 


1.02 (froz.) 


2.82 (froz.) 


4.3±1.6 


2.2±0.4 


0.21 (froz.) 


0.88 


0.72 



* Chandra data (total exposure time) 

Table 3. Best-fit spectral parameters for V892 Tau and V892 Tau NE during the Chandra observation and the second 
XMM-Newton exposure. The spectral parameters are derived assuming a single-temperature optically thin plasma 
model. 



Source 


jV(H) 


kT 


EM 


Z 


x 2 


P 




10 22 cm" 2 


keV 


10 oa cm" d 


Z 






V892 Tau (Chandra) 


0.83±0.08 


2.10±0.19 


4.31±0.95 


0.06±0.07 


1.04 


0.40 


V892 Tau (XMM f) 


0.65±0.06 


1.53±0.16 


5.19±1.49 


0.09±0.05 


1.33 


0.06 


V892 Tau (XMM 2) 


0.97±0.05 


8.11±1.00 


21.22±1.01 


0.16±0.09 


0.90 


0.79 


V892 Tau (XMM 3) 


1.00±0.04 


6.68±0.58 


26.44±1.11 


0.26±0.08 


1.01 


0.46 


V892 Tau NE 


1.10±0.21 


1.08±0.22 


0.78±0.56 


0.2 (froz.) 


0.40 


0.98 



XMM 1: X-ray quiescent, XMM 2: rising phase of large flare, XMM 3: ffare maximum 



As already mentioned, V892 Tau and V892 Tau NE are 
unresolved by the XMM-Newton PSF, nevertheless the ev- 
idence that the origin of the observed flare is the Herbig 
Ae star is compelling. First, as summarised in Tabled the 
XMM-Newton position, before and during the flare, of the 
source associated with system the V892 Tau is between 
2.4 and 3.f arcsec from the radio position of the main star 
V892 Tau, while is between 6 and 6.8 arcsec off the posi- 
tion of V892 Tau NE. The absolute measurement accuracy 
of the XMM-JVewton pointing is 4 arcsec (XMM-JVeraton 
Users' Handbook), so the offset of the source in the EPIC- 
pn image from the radio coordinates of V892 Tau is well 
within f er of the pointing accuracy. This is fully consistent 
with V892 Tau being the dominant source of X-ray emis- 
sion throughout the XMM-Newton observation, as it is the 
case during the Chandra observation, where V892 Tau is 
roughly 10 times more luminous than V892 Tau NE. 

Second, the XMM-Newton position of the source asso- 
ciated with the system V892 Tau before and during the 
flare is constant within 0.7 arcsec , while the angular sep- 
aration of V892 Tau and V892 Tau NE is 4.1 arcsec. At 
the same time (as shown in the previous section) the spec- 

4 The 0.7 arcsec shift of position of the source before and 
during the flare is to the south, while V892 Tau NE is to the 
north-north east of V892 Tau. Therefore, although not signif- 
icant given the ~ f arcsec accuracy of XMM-Newton relative 
pointing error, the position shift is consistent with V892 Tau 
becoming even stronger in respect with V892 Tau NE. 



tral characteristic of the EPIC-pn data on the V892 Tau 
system while quiescent are in good agreement with the 
characteristics derived from the spectral data on V892 Tau 
from Chandra but significantly different from the ones de- 
rived from the Chandra spectral data on V892 Tau NE. 
These two facts strongly indicate that the XMM-Newton 
data are, for both the quiescent emission and the flare, 
fully dominated by the emission of V892 Tau. Had the 
flare been associated with V892 Tau NE, a shift of the 
X-ray source position during the flare would have been 
observed. 

Finally, if one wanted to explain the XMM- Newton 
light curves as mostly determined by V892 Tau NE, dis- 
regarding the evidence coming from the pointing and 
spectral information, one would have to require that the 
source which does not appear significantly variable in the 
Chandra observation becomes from ten (before the flare) 
to a hundred times (during the flare) brighter during 
XMM-Newton observation. This is unlikely. 

Given the facts above, we conclude that the XMM- 
Newton light curve shown in Fig. 21 an d the EPIC-pn 
spectra shown in Fig. El and Fig. are determined by the 
activity of the Herbig Ae star V892 Tau and minimally in- 
fluenced by the presence of V892 Tau NE. Consequently, 
below we discuss the XMM-Newton data as originating 
from V892 Tau and do not further comment on the unre- 
solved presence of V892 Tau NE. 
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Fig. 6. Background-corrected spectra of the V892 Tau 
system during the tree different phases of the second 
XMM- Newton observation. From top to bottom: 14 ks 
integration time before the flare, 8 ks during the rising 
phase of the flare, 10 ks during the flare maximum. 



4.1. Quiescent emission 

By "quiescent" here we mean the time when V892 Tau is 
not undergoing the large flare, even though during this 
phase we do observe smaller flare-like events (both in 
Chandra and XMM- Newton data). 

The Chandra spectrum and the XMM-iVewton spectra 
of V892 Tau while the source is quiescent can all be well 
described by an absorbed two temperature plasma model 
with N{B) ~ 0.9 x 10 22 cm" 2 kT x = 1.0 keV, kT 2 = 



z.8 kev and a relative metal abundance oi Z = O.z Zq. 
As discussed above, for the Chandra data a slightly higher 
value of iV(H) is required in order to fit the spectrum, con- 
sistent with the known contamination of the ACIS chips. 
Here therefore we only refer to the value for iV(H) derived 
from the XMM-Newton data. 

The X-ray emission from V892 T au has been 
studied previously by IStrom fe Stroml l|l994h and 
IZinnecker fc Preibischl ll 19941) . These two works use the 
same ROSAT obs ervation but derive fr om it different 
model parameters. IStrom fc Stroml l|l994f) fit the spectral 
data with an absorbed two temperature plasma model 
with 2V(H) ~ 1.3±0.6x 10 22 cm- 2 , kT x = 0.55, kT 2 = 1.20 
and estimate an X-ray lumi nosity Lx = 1-0 x 10 31 e rg s 



(in the 0.2-2.4 keV band). IZinnecker fc Preibischl (Hi 
model the spectra of V892 Tau with an absorbed 1 temper- 
ature plasma with N(R) = (4.8 ± 1.1) x 10 21 cm" 2 , kT = 
2.3±0.4keV and estimate L x = (1.9 ± 0. 8) x 10 30 ergs" 1 
(a fac tor of ~ 5 lower than the estimate of lStrom fc Stroml 

We derive iV(H) ~ 0.9 x 10 22 cm " 2 , compatible to 
the estimate of lStrom fc Stroml ||1 9941) . This value for the 
absorbing column density corresponds to a visual extinc- 
tion 5 of Ay ~ 4.7, which is cl ose to the value estimated by 
IZinnecker fc Preibischl ll 19941) from the B — V value given 
in iHerbig fc Bell i fa9§8l) an d a factor of ~ 2 lowe r than 
derived hv lFliasTl|l97Sl) and IStrom fc Stroml lll994l) 6 . We 
note that a lower value of visual extinction is consistent 
with the star spectral type being A6 rather than B9. 

The two temperatures for the two component plasma 
model that we der ive here are significan tly higher than 
the values given bv lStrom fc Stroml (|l994h : this is not un- 
expected given the much softer (and nar rower) bandpass 
of the ROSAT PSPC instrument used bv lStrom fc Stroml 
l|l994|) . The temperature derived from our 1 plasma com- 
ponent fit to the Chandra data (kT = 2.1 keV) is how- 
ever consistent with plasma temperature estimated by 
IZinnecker fc Preibischl l|l994h . The values we derived for 
the two plasma temperatures in V892 Tau are somewhat 
higher than the typical values derived for the X-ray emis- 
sion from low mass PMS in Taurus: in an analysis of the 
spectral characteristic of 9 T Tauri sta rs in L1551 we foun d 
typically kT x - 0.3 and kT 2 ~ 1.2 (jFavata et al.ll2003h . 
On the other hand the value for the plasma metallicity 
Z = 0.2 for V892 Tau is typical of the value we derived in 
the same study for the 4 Weak Lined T Tauri stars in the 
sample. 

The intrinsic luminosity we estimate for V892 Tau 
in its quiescent state is Lx = 1.6 x 10 30 erg 
s" 1 , in good agreement with the value given by 
IZinnecker fc Preibischl l|l994|) and a factor of ~ 6 lower 



using 
atoms cm 



conversion 
-l 



1.9 x 10 2 



mag- 



factor mm/A v 
see e .g. lCoJ200nl) . 
e IStrom fc Stroml dl994) estimate Ay from optical observa- 
tion. This may explain why their value for the absorbing col- 
umn density is close to our estimate while their estimate for 
Ay is a factor of 2 higher. 



than derived by Strom & btrom ( 1994). Ihis indicate that 
IZinnecker fc PreibischTijl994|) are probably correct when 
they suggest that the reason for t he discrepancy between 
their values and the one derived bv lStrom fc Stroinl (^994) 
is the fact that Strom & Strom may have erroneously in- 
cluded the nearby source [BHS98] MHO 11 in the source 
circle of V892 Tau. 

A wind-related origin of the X -ray emission has been 
propo sed f or Herbig Ae/Be stars bv lzinnecker fc Preibischl 
l)l994[) and lDamiani et all l)l994[) . This scenario, however, 
seems an unlikely explanation for the origin of the X-ray 
emission on V892 Tau in its quiescent phase. 

As described in Sect. during this phase, V892 Tau 
presents significant short term variability, with its X-ray 
flux varying by a factor of 2 in less than 1 ks. The impulsive 
rise in flux appears to be associated with an hardening of 
the spectrum in the Chandra data. This type of variabil- 
ity is similar to the one observed in lower-mass pre-main- 
sequence stars, where the X-ray emission is of coronal ori- 
gin, while being substantially different fro m the 20-30% 
X-ray flux variations observ ed in OB stars ijCollura et alJ 
Il989t ICassinelli et aliHool where the emission mecha- 
nism is wind-related. 

In addition, from the l uminosity of V892 Tau esti- 
mated to be around 38 Lq ijBerrilli et al .111 9921) . its ratio 
of X-ray luminosity to the total luminosity while quiescent 
is Lx/^boi — 1 x 10~ 5 . This is two orders of magnitudes 
greater than the typical ratios found for OB stars. On the 
other hand a value of Lx/iboi — 10~ 5 is not uncommon 
for low mass stars in which the X-ray emission mechanism 
is coronal. 

The possibility of a wind related origin of the X-ray 
emission from V892 Tau in its quiescent phase appears 
therefore unlikely. This is in agreement with the conclu- 
sions reached in the two statistical studies of the prope rties 
of Herbig Ae/Be s t ars of|Preibisch fc Zinneckerl (|lp96) and 
lHamaguchi et al 1 ll2002aF These studies indicate that the 
X-ray emission from Herbig Ae/Be stars is generally asso- 
ciated with higher plasma temperature and higher X-ray 
to bolomctric luminosity ratios than typically observed for 
the X-ray emission from main sequence OB stars. 

A wind scenario, finally, cannot account for the ob- 
served flare event of V892 Tau. 

4.2. Flare event 

During the flare event the X-ray luminosity of V892 Tau 
increases by a factor of ~ 15, from Lx = 1.6 x 10 30 erg 
s _1 to Lx = 2.4 x 10 31 erg s" 1 , while the temperature of 
the plasma increases from kT = 1.5 keV to kT — 8.1 keV. 
The source luminosity increases over a relatively long time 
of ~ 10 ks, and hovers close to the luminosity maximum 
for at least another 10 ks thereafter - the end of the ob- 
servation does not unfortuna tely allow to study the decay 
phase. As demonstrated by iReale et"aT. | ll2002l) flares of 
this type cannot take place in the freely expanding plas- 
moids of a stellar wind. The observed slow increase of X- 



ray luminosity at the beginning ot the tiare event requires 
the presence of a confining magnetic field. 

In order to gain some quantitative insight on the event 
and given the similarities of the derived light curve with 
the ones of other stellar flares we have analyzed the 
event through scaling obtained from detailed hydrody- 
namic models of flaring plasma confined in a closed coronal 
loops, as in solar flares. 

4.2.1. Flaring region characteristics 

It is customary to derive information about the size of the 
flaring loops from the flare evolution, i.e. the light curve. 
It has been shown that the decay time of the light curve is 
linked to the plasma cooling times, which, in turn, depends 
on the leng th of the loop which confines the plasma (e.g. 
lReak-jl2f)02l and references therein): the slower the decay, 
the longer the loop, unless a significant residual heating 
sustains the decay and makes the decay-time/cooling-time 
dependence less tight. 

The XMM-Newton observation of the flare on V892 
Tau does not cover the decay phase and therefore diag- 
nostics using the characteristics decay times are not feasi- 
ble. For this particular flare, however, we are able to infer 
some information on the size of the flaring region from the 
observed rise phase. 

The evolution of the flaring plasma confined in a loop 
is well-k nown from extens ive hydrodynamic loop model- 
ing (e.g. lPeres et al.lll9 82): a strong heating pulse is trig- 
gered in an initially quiescent coronal loop and makes the 
temperature increase by up to several tens of MK along 
the whole loop in a few seconds, due to the high plasma 
thermal conduction. The dense chromosphere at the loop 
footpoints is heated violently and expands upwards with 
a strong evaporation front. The rising plasma fills up the 
loop, very dynamically first and then more gradually, ap- 
proaching a new hydrostatic equilibrium at a very high 
pressure. The loop X-ray emission increases mostly fol- 
lowing the increase of emission measure, and forms the 
rise phase of the flare. 

Although the evolution in the rise phase is very dy- 
namic and non-linear, we can nevertheless derive an ap- 
proximate time scaling. From the equation of energy con- 
serva tion of the confined plasma (e.g. Eq.(3) in lSerio et alJ 
19911) it can be seen that, after the very initial seconds, 
dominated by the plasma kinetics, the evolution in the 
bulk of the rise phase can be approximately described as a 
linear increase of the plasma internal energy density driven 
by the (constant) energy input rate per unit volume: 



The pressure increases linearly as well, with a rate 
p' ~ 2/3Eh- The linear trend goes on until other ener- 
getic terms, and in particular the plasma radiative losses, 
become important in the energy balance. Then the pres- 
sure increase slows down to approach asymptotically the 
equilibrium condition described by the loop scaling laws 



IjKosner et al.H197q) tor a loop with heating volume rate 
Eh- A constantly linear pressure increase with a rate p' 
would imply that the loop reaches the equilibrium pres- 
sure in a time t p w po/p' = 3/2po/Eh- 

Detailed hydrodynamic simulations of the flaring loop 
rising phase have been performed on purpose for the 
present work. We simulated the evolution of the plasma 
confined in a coronal loop brought to flare conditions by 
a strong impulsive heating, as done in previous extensive 
modeling of solar and stellar coronal flare (Peres et al. 
1982, 1987; Reale et al. 1988; Reale & Peres 1995; Betta 
et al. 1997, 2001). The assumption specific for this work 
is that of a constant flare heating lasting long enough 
to let the loop reach equilibrium conditions. The simula- 
tions show that the pressure linear increase rate is ~ 20% 
slower than p' . If we also include the slowing of the pres- 
sure increase at later times, we obtain that the time taken 
to reach ~ 90% of the equilibr i um pr essure is twice t p , 
i.e. dtp w 3pq/Eh- ISerio et all l)l99l|) have shown that 
T s = Po/Eh ~ 120ig/v / T7 is the loop thermodynamic 
decay time, where Lg is the loop half-length in units of 
10 9 cm and Tj is the plasma maximum temperature in 
units of 10 7 K. Therefore, we end up with an expression 
St p ~ 3rs which links the time taken by the loop to ap- 
proach the flare equilibrium condition to the loop thermo- 
dynamic decay time. 

We note that in this flare the rise phase is relatively 
long (~ 8 ks) , and the flaring loop may therefore approach 
equilibrium condition. So, if we set St p equal to the rise 
duration time of ~ 8 x 10 3 s, we obtain an estimate for 
Lg, once T7 is known. In this case T 7 - 18 7 , so that we 
obtain Lg ~ 100. 

This estimate for the size of the flaring region has been 
obtained under a few assumptions which deserve some 
comments. First, we have considered a flare occurring in 
a single loop. This may not hold true for such an intense 
and long-lasting flare, which may perhaps be described 
more properly as a two-ribbon flare, consisting of progres- 
sively reconnecting higher and higher loops. Nevertheless, 
it often occurs in solar two-ribbon flares that the rise phase 
mostly involves a dominant loo p structure, and then grad- 
ually extends to others (e.g. lAschwanden fc Alexander! 

We may associate the estimated length to such 
dominant structure. 

Another implicit non-trivial assumption is that of 
a heating pulse constantly high during the rise phase. 
Indeed, a gradually increasing heating function may drive 
the observed gradual rise of the light curve, invalidating 
the estimations made above. However, in such a case we 
should observe also a gradual increase of the temperature, 
while the indications are for a sudden jump of the temper- 
ature to the flare value, which is more typical of a heating 



7 iBriggs fc Pvd ll200.?l) give the empirical expression T max = 
0.184 Jobs 3 linking the temperature obtained from spectral fit- 
ting of EPIC data (T b 8 ) to the loop maximum temperature 
(T m 

ax ) 



pulse. We may be therefore quite confident that the total 
loop length is ~ 2 x 10 11 cm. 

If we assume a loop aspect R/L ~ 0.1, where R is 
the radius of the loop cross-section, assumed circular and 
constant along the loop, we obtain a total loop volume 
V ~ 6 x 10 31 cm 3 ; a maximum emission measure of ~ 2.6 x 
10 54 cm -3 then implies a maximum average loop plasma 
density of ~ 2 x 10 11 cm -3 and a maximum pressure of 
the order of 10 dyne cm -2 . This value is compatible with 
the equilibrium pre ssure obtained from loop scaling laws 
llRosner et al.lll978|) and consistent with the hypothesis 
of the loop at maximum X-ray luminosity being close to 
equilibrium conditions. In order to confine a plasma at 
such a pressure, a magnetic field of more than ~ 500 Gauss 
is required. The origin of such a field in stars that are 
thought to be fully or nearly-fully radiative is puzzling. 

4.2.2. What is the origin of the magnetic field? 

In low-mass stars the presence of a significant convection 
zone supports the dynamo mechanism that can generate 
the confining magnetic fields at the origin of the X-ray 
flare events, but according to classical models, pre-main- 
sequence star with masses in excess of 2 Mq are expected 
to follow fully radiative tracks once the quasi-static con- 
traction has ended. 

IPalla fc Stahlerl (1990) have made the suggestion that 
the surface activity and winds observed in Harbig Ae/Be 
may be related to the presence of an outer convection zone. 
In their interpretation, this convection zone results from 
the subsurface shell burning of residual deuterium which 
was accreted d u ring the protostar phase. Nevertheless 
IPalla fc Stahlel |l993) reconsider their hypothesis, and 
explain that according to their models the retreat of the 
proto-star outer convection zone does last a substantial 
fraction of the pre-main-sequence lifetime of an interme- 
diate mass star. During this retreat, however, the effective 
temperature remains relatively low, so that the star would 
not be observed with an A or B spectral type. They con- 
clude that the presence in Herbig Ae/Be stars of surface 
activity and strong winds is not linked to an outer con- 
vection zone, since their model shows that such convection 
always vanishes with the rising effective temperature. 

Recently ISiess et all l)2000j) have presented calcula- 
tions of pre-main-sequence evolutionary tracks for low- 
and intermediate-mass stars. These models predict the ex- 
istence of a thi n convective enve l ope in young AB stars. 
In their review, iFavata fc Micelal l|2003|) suggest that this 
thin convective envelope, of roughly 2 x 10~ 3 times the 
stellar radius, could be at the origin of a low-coronal ac- 
tivity (at the level of the observed minimum for solar type 
stars) in Altair (A7V) and thus explain the source X-ray 
luminosity of Lx = 3 x 10 27 erg s _1 . 

According to the same models a star with mass 
1.9 Mq, the same metallicity of the Sun and an age of 
10 million years, would have a spectral type A6, a lumi- 
nosity of 12 Lq, not too far from the one of V892 Tau, 



and a convective envelope oi a traction ot 0.0018 its radius. 
Since the model predicts a stellar radius of 1.2 x 10 11 cm, 
the size of the convective region would be 2.2 x 10 8 cm 
(~ 0.3% Rq). It is not clear whether such a thin con- 
vective envelope, that may be sufficient to generate the 
low coronal activity invoked to explain the 3 to 4 order 
of magnitude fainter X-ray emission of a main sequence A 
star, can sustain the dynamo action necessary to explain 
the strong X-ray activity of V892 Tau. Indeed, from our 
flare model we derive a flare loop length of ~ 2 x 10 11 cm. 
This is comparable to the stellar radius and corresponds 
to a size of ~ 500 times the thin convective envelope. 

An alternative mechanism to sustain the dynamo ac- 
tivity in these predominantly ra diative stars has been pro- 
posed bv lTout fc Pringlel l|l995^ They argue that dynamo 
activity can be sustained in AB stars for a substantial 
fraction of their pre-main-sequence life time by tapping 
the initial stellar differential rotation - or shear energy. 

5. Conclusions 

We have analysed the light curves and spectral data of the 
system V892 Tau and V892 Tau NE in a Chandra 18 ks 
exposure and 2 consecutive XMM-Newton exposures of 74 
and 45 ks (nominal) . In the Chandra data, the Herbig Ae 
star V892 Tau is well resolved from the low mass later type 
apparent companion V892 Tau NE. During the Chandra 
observations V892 Tau shows significant variability, with 
its X-ray flux varying by a factor of 2 in less than 1 ks. 
This type of variability is reminiscent of the type of flaring 
variability observed in lower-mass pre-main-sequence star 
and has never been reported before for a Herbig Ae star. 

During the second XMM-Newton exposure of the sys- 
tem V892 Tau and V892 Tau NE a large flare event takes 
place. The source luminosity impulsively increases by a 
factor of 15, from 1.6 x 10 30 erg s _1 to 2.4 x 10 31 erg 
s , while the temperature of the plasma increases from 
kT = 1.5 keV to kT = 8.1 keV. 

V892 Tau and V892 Tau NE are unresolved by the 
XMM-Newton PSF, nevertheless the combined XMM- 
Newton and Chandra data set provide strong evidence 
that the origin of the observed flare is the Herbig Ae star. 

We have modeled the flare event and find that a mag- 
netic field of 500 Gauss in intensity is necessary to con- 
fine the plasma which reaches a temperature in excess of 
100 MK. The generally accepted mechanism that can sus- 
tain confining magnetic fields of this intensity is a dynamo 
action. Therefore under the assumption that surface con- 
vection is a necessary ingredient for a dynamo, our findings 
imply the presence of a convective envelope in the Herbig 
Ae star V892 Tau. 

Recent models by ISiess et all l)200fll indicate that a 
pre-main-sequence A6 star similar to V892 Tau has a very 
thin convection zone, extending to a depth of only 0.2% of 
the stellar radius. Whether this thin convective envelope is 
sufficient to sustain the dynamo action required to explain 
the vigorous X-ray activity of this Herbig Ae star remains 
to be assessed. 
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